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I. INTRODUCTION 

Gravitational waves (GWs) astronomy opens a new 
window to study the physical processes in the very early 
universe: relic GWs propagate almost freely throughout 
the universe expansion, and thus they retain the infor¬ 
mation about the physical conditions and physical pro¬ 
cesses at the moment of their generation (see for reviews, 
[in and references therein). There are various mech¬ 
anisms that might generate such GWs. In the present 
paper we focus on the generation of GWs during cos¬ 
mological electroweak (EW) and Quantum GhromoDy- 
namic (QGD) phase transitions (PTs) through the tur¬ 
bulent helical sources which can arise and follow the PT 
bubble collisions. The GW generation mechanism asso¬ 
ciated with bubble collisions during the first order PTs 
has been widely discussed in literature, starting^from the 
pioneering works ill and re-addressed later [3-1^ 

For a cosmological phase transition to produce strong 
enough turbulent motions and magnetic fields, which will 
result in the detectable signal of GWs, they must be first 
order PTs, with bubble formation and bubble collisions. 
For the EWPT, with the standard EW Lagrangian plus a 
Stop, the supersymmetric partner of the top quark, called 
the MSSM EW Lagrangian, the EWPT is first order PT 
[ 2 ^. The QCDPT has been shown to be first order PT 
by lattice gauge calculations [13, . 

Both turbulent motions and matic fields can produce 
relic GWs through their anisotropic stresses, see Refs. 
(29l - l39j| . It has been pointed that the GWs generated 
by magnetic fields can be detected through Laser Inter¬ 
ferometer Space Antenna (LISA) [40l - l43| . In difference 
to the GWs sourced solely by PTs bubble collisions, the 
presence of turbulent (kinetic and magnetic sources) in¬ 
creases the detection pr osp ects [dsj not only from EWPT 
but from QCDPT too [d^.^ One of main goals of Euro¬ 
pean Space Agency (ESA) - NASA planned join mission 
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^ GWs from QCDPT are potentially detectable through pulsar 
timing, see Ref. [ 4 ^ and references therein. 


LISA [ 4 ^, was the detection of low frequency GWs (sub- 
Hz region). The new development of this program is the 
European only ESA mission, so called New Gravitational 
wave Observatory (NGO) - aka eLISA (evolved LISA) 
( 4 ^ . One of major parts of its science program consists 
on the direct detection of GWs from cosmological PTs, 
see Refs, [sol - ls^l for details. 

In the present paper we extend our previous study Ref. 
[s^ , and we investigate the degree of polarization of GWs 
generated via cosmological PTs through helical hydro 
and magnetized turbulent sources using the formalism 
given in Ref. [s^. We adjust the previous formalism to 
determine the polarization degree of GWs from helical 
kinetic turbulence to a more complex scenario of MHD 
turbulence present during the cosmological PTs. More 
precisely, we use the recent results of numerical simula¬ 
tions [54 - 1^ to set the statistical properties of helical 
MHD turbulence. Another difference from the formalism 
of Ref. [s^ consists in computing the energy density and 
peak frequency of GWs using the analogy with acous¬ 
tic waves production by hydrodynamical turbulence [s^ 
(which we can call aeroacoustic approach [53 - l59l| ). 

Gharge-conjugation-Parity (CP) violation is necessary 
for the production of magnetic helicity via bubble col¬ 
lisions, [ 13 . EWPT and QCDPT bubble collisions re¬ 
sult in development of helical (kinetic or/and magnetic) 
turbulence, due in part to CP violation, which will lead 
to circularly polarized GWs back grou nd. In the case of 
strong enough helical sources [^ l.'Ifill . the degree of po¬ 
larization is potentially detectable |6ll - l^ .^ 

In our present study we follow the helical (chiral) mag¬ 
netic fields generation scenarios (during PTs through 
bubble collisions) presented in Refs. [ll, (EWPT) 
and Refs. Hill (QCDPT) (see also Ref. [H for a brief 
review of these models). Upon generation the magnetic 
field starts to interact with primordial plasma that leads 
to development of magnetically dominant MHD and sec- 


^ The indirect tool to detect circularly polarized GWs consists 
on searching parity violating signals on cosmic microwave back¬ 
ground maps, see Refs [6347111 for original studies and Ref. 
for a review. 
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ondary kinetic turbulence, for pioneering studies see Refs. 
(TGI - IT^ . In what follow we adopt the results of numerical 
simulations of Refs. ^ and their phenomenologi¬ 

cal interpretation given in Refs. Bin 

The structure of the paper is as follows: In Sec. II we 
review the GWs generation formalism and define the cir¬ 
cular polarization degree of GWs. We discuss the hydro 
and MHD helical turbulence modeling in Sec. Ill and 
compute the GW signal and its polarization in Sec. IV. 
We give our results for both EWPT and QCDPT gener¬ 
ated GWs in Sec. IV, and we conclude in Sec. V. We use 
natural {h = \ = c) Lorentz-Heaviside units. 


II. GRAVITATIONAL WAVES GENERATION 
OVERVIEW 


where hij (x, t) is the tensor metric perturbation, t is 
physical time, i and j are spatial indices (repeated in¬ 
dices are summed), and G is the gravitational constant. 
We have neglected the term oc dhij{x,t)/dt due to our 

(T) 

assumption of the short duration of turbulence. Pb 
{the script ”T” indicates that we are interested on the 
tensor part of the turbulent source) is the traceless part 
of the stress-energy tensor (x,t), which is constructed 
from kinetic (K) or magnetic (M) turbulence normal¬ 
ized vector fields^ (as we will show below the equipar- 
tition is established between kinetic and magnetic turbu¬ 
lent motions which simply doubles the source term, i.e. 


n: 


(K) 


n 


(M) 


~ 2 n[p) given by [si) : 


(x, t) = Ty (x, t) - i dy r(x, t), (4) 


We assume that GWs are generated through kinetic 
and MHD turbulence which follow the PT bubble colli¬ 
sions [3, l29l - l^ [s^. To be general as possible we present 
the common description for EWPT and QCDPT, dehn- 
ing the PT temperature as T* (T* = 100 GeV for EWPT 
and T* = 0.15 Gev for QCDPT, and the typical proper 
length scale through Iq which can be associated with the 
PT bubble size lb (the assumption Iq ~ lb is well justi¬ 
fied because in our theory bubble collisions during PT 
generate a magnetic field at bubble walls, and this ini¬ 
tial field starts to interact with primordial plasma result¬ 
ing in development of kinetic and MHD turbulence with 
a typical length scale that corresponds to the magnetic 
field injection scale) [s^. We also define 5 * as a num¬ 
ber of relativistic degrees of freedom: for the standard 
model we have g* = 106.75 as T —>• 00 . (g* = 100 for 
EWPT and g* = 15 for QCDPT). In our further con¬ 
sideration we assume that the duration of the turbulent 
sources, tt are short compared to the universe expansion 
time-scale at PTs, i.e. tt < with H~^ the Hubble 
radius at PT. This assumption makes possible to neglect 
the expansion of the universe, although limits our con¬ 
sideration by the GWs background generation only from 
PT, and completely neglects GWs arising from decaying 
turbulence (which might last log-enough after the end of 
PTs). 

GWs (the tensor metric perturbations above the 
standard Friedmann-Lemaitre-Robertson-Walker homo¬ 
geneous and isotropic background) are generated from 
turbulence (including kinetic and magnetic fluctuations) 
through the presence of anisotropic stresses as 

^2 

V^/iy (x,t) - —/iy (x,t) = -167rGn|J^(x,t), (1) 


^ We underline the nature of the secondary character of fluid mo¬ 
tions, because the bubble collision itself might lead to develop¬ 
ment of purely hydrodynamical turbulence during PTs, see Refs. 

, while here we note that bubble collisions result in generation 
of magnetic fields 


here T = [T]^ is the trace of the Tij tensor. 

As we can expect the kinetic and magnetic turbulent 
fluctuations generate stochastic GWs, which can be char¬ 
acterized by the wave number-space two-point function 
as, 


{h*j{k,t)him{'k',t + t)) = ( 27 r)^d(^^(k - k') x 


X 


Mijim{k)H{k, t) + iAijim{k)'H{k, r) 


.(5) 


Here we use the Fourier transform pair of the ten¬ 
sor perturbation as: /iij(k, t) = J d^x e^^'^hij(xA) and 
= f d^ke~^^'^hij(k,t)/(27r)^. The brackets 
(...) denote an ensemble average over realization of 
the stochastic source. The spectral functions H{k,t) 
and 'H(fc,t) determine the GW amplitude and polar¬ 
ization, 4Mijlm{k) = Pil{k)Pjm{k) + Pim{k)Pji{k) - 
^ij {^^^lm{k^ J and SAijlm{k') = kq^Pjjfi{k{cilq 
Pii{k)ejmq + Pi7n{k)ejiq -f Pji{k)eimq] are tensors, with 
the projection tensor Pij{k) = Sij — kikj (with 6ij - the 
Kronecker delta, kt = ki/k and k = |k|), is is the 
totally antisymmetric symbol. We choose GW prop¬ 
agation direction pointing the unit vector 63 , and we 
use the usual circular polarization basis tensors efj = 
— (ei ± ie 2 )i X (ei ± ie2)jlV2. We define two states 


^ The kinetic and magnetic perturbation stress-energy tensor are 

=wui(x,t)uj(x,t), (2) 

= wbi{x,t)bj{x,t), (3) 

where w = p -h p is enthalpy of the fluid with density energy, p, 
and pressure p, u(x, t) is the kinetic motion velocity field and 
b(x, t) is normalized magnetic field, b = B/V^ttw, that repre¬ 
sents the Alfven velocity, va of the magnetic field. The normal¬ 
ized energy of the magnetic held is then £m{v) — 
while the normalized kinetic energy density is given through 
SK{t) = {u^(£))/2. The advantage of such representation con¬ 
sists on eliminating the expansion of the universe, since physical 
and comoving values of the normalized magnetic field amplitude 
are the same. 
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and h~ corresponding to right- and left-handed circularly 
polarized GWs, = h+e^- -I- h~e~y Through above no¬ 
tations the circular polarization de gree was derived for 
GWs from Gamma Ray Bursts by [8^ and in the con¬ 
text of cosmological GWs is reproduced [s^, 

GW.- ^ {h+*{k)h+(k')-h-*ik,)h-ik')) ^ 

^ ^ {h+*{k)h+{k') + h-*{k)h-{k')) H{k)' 

( 6 ) 

Here we omit the time dependence of the polarization 
degree 

As we already underlined we are interested on GWs 
generation only from a short duration sources acting dur¬ 
ing PTs. After generation the GWs propagate almost 
freely, and we account for the expansion of the universe 
by a simple re-scaling of the frequency and the amplitude 
by a factor equal to 

^ . 8 X 10- f 

Oo V 



This factor is safely canceled when computing 
although the more complex consideration of decaying tur¬ 
bulence (long lasting sources) will make {k) time 
dependent function. 

To estimate the polarization degree of GWs from PT 
generated helical fields we need to compute two spectral 
functions 'H{k) and H{k) at the moment of PT, which 

are determined by the helical anisotropic sources 

and . In our previous work we have computed the 
typical amplitude and frequency of GWs generated dur¬ 
ing cosmological PTs [s^. In Ref. [s^ the polarization 
degree of GWs from kinetic (hydro) turbulence has been 
estimated. It has been shown that fully helical turbu¬ 
lence leads to {k) —>■ I. In the present work we 
follow the GWs generation formalism from helical mag¬ 
netized sources presented in Refs. and apply it 

to the EWPT and QCDPT cases. 


III. KINETIC AND MHD TURBULENCE 
MODELING 


The magnetic field amplitude (i.e. total magnetic field 
energy density) is strongly limited by the big bang nucle¬ 
osynthesis (BBN) bound requesting that the total mag¬ 
netic field energy density cannot exceed the 10% of the 
radiation density at the moment of the magnetic field 
generation. In terms of the effective comoving magnetic 
field value, Reff — 8.4- 10“'^(100/(7*)^/® Gauss (G) or in 
the terms of Alfven velocity va = |b| < 0.4 [s^ . 

As we noted above the magnetic field generated at 
one scale (the magnetic field initial spectrum can be 
approximated being peaked at the typical wavenumber 
fco = 27r//o, so in Fourier k-space described by (k— kp) 
function) after interactions with plasma leads to devel¬ 


opment of turbulence,® and the sharply peaked initial 
spectrum is redistributed respectively. For isotropic sta¬ 
tionary turbulence the normalized magnetic vector field 
two-point correlation function is: 

{b*{k)b,{\f)) = (2^)3d(3)(k - k') F^{k) (8) 

where 

F^(k) = P,j{k)SM{k) + ie,,ikiAM{k). (9) 

The power-law spectral function SMik) = Sok'^^ and 
Anik) = determine the energy density 

and current helicity of the magnetic field® and ns and 
UA are magnetic field and helicity spectral indices which 
determine the spatial distribution of the magnetic field 
and its helicity. The establishment of stationary turbu¬ 
lence with the stationary (time- independent) two-point 
correlation function given through Fqs. ® ® requires 

the presence of long-lasting sources. To account for the 
short acting PT turbulent source (the turbulence dura¬ 
tion time tt is short enough compared to the universe 
expansion time-scale we have to modify the spec¬ 

tra Smik) and AM{k) making them time dependent, see 
below. 

Following the description of hydro- and MHD turbu¬ 
lence generated during PTs, we distinguish three spa¬ 
tial spectral sub-regimes of turbulent fluctuations: (i) 
the large scale decay range kn^ < k < ko (where phys¬ 
ical wavenumbers kn^, = 27r/i7~^ and ko correspond to 
the PT Hubble length scale and the largest PT length 
size); the minimal wavenumber corresponds to the Hub¬ 
ble scale beyond which causally generated magnetic 
field is frozen-in and any interactions are forbidden due 
to causality requirement; (ii) the turbulent (or so called 
inertial) range ko < k < kjj (where kjj the damping 
scale of turbulence through viscous dissipation and mag¬ 
netic resistivity, which is determined by plasma proper¬ 
ties); (iii) the damping range k > ku. All these typical 
wavenumbers ko, and ko) are time-dependent due 

to interactions of magnetic field with plasma and the ex¬ 
pansion of the universe. Magnetic helicity presence plays 
here a crucial role leading to re-arrangement of the helical 
structure at large scales [s^- The expansion of the uni¬ 
verse lead to additional effects: namely, the PT bubble 
size determined length scale lo is strengthened by a fac¬ 
tor a{t)/ai, (being oc during the radiation-dominated 


® Primordial plasma is a perfect conductor with extremely high 
values of kinetic and magnetic Reynolds numbers, and current 
numerical simulations are still behind to approach necessary reso¬ 
lutions and timescales to describe adequately physical conditions 
and processes in the early universe. 

® Magnetic helicity defined as {A(x) ■ B(x)) is a gauge-dependent 
quantity, while normalized (or regular expressed through B(x)) 
current helicity {b(x) ■ [V X b(x)]) is gauge independent, see Ret. 
[s^ for details, also it allows the direct analogy with the kinetic 
helicity {u(x) ■ [V X u(x)]), and thus to consider both helical 
sources in a common formalism. 
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epoch and oc during the matter-dominated epoch), 
while the Hubble length scale oc t. As a result a per¬ 
turbation with kHo < k < kn^, (with kno = 2tt/Hq and 
Hq is the today Hubble radius) will enter the Horizon at 
some point [83 • Since we are focused on the short dura¬ 
tion sources, we will completely neglect the GWs signal 
for the large-scale decay range k < ko (where the spectral 
shape of the field is given through the causal Batchelor 
spectrum with ns = 2 [13). Obviously the GWs signal 
from the viscous damping range k > kjj is also negligibly 
small. 

The realizability condition implies that |AM(fc)| < 
SM{k) (the modulus sign reflects a possibility of having 
positive or negative helicity). The spectral indices values 
ns and ua strongly depend on the turbulence model. 
In the inertial range (/cq < k < ku), for non-helical 
turbulence the Kolmogoroff model implies ns = —11/3. 
Some models lead to the different spectral shapes such as 
ns = —7/2 - Iroshnikov-Kraichnan model for magnetized 
turbulence ns = —4:- the weak turbulence model 

or magnetically dominant turbulence [o^. In the 
presence of helicity the consideration is even more com¬ 
plex, and requires careful investigation through numeri¬ 
cal simulations that it is beyond the scope of the present 
paper. Based on the phenomenological dimensionless de¬ 
scription, if the process is driven by the magnetic energy 
dissipation at small scales, it is assumed ns = —11/3 and 
nA = —14/3 (so called helical Kolmogoroff model) [^ . 
while if the process is determined by helicity transfer (in¬ 
verse cascade) and helicity dissipation at small scales it 
is adopted ns = —13/3 = n^ [^ . 

To account the short-duration turbulence (not enough 
to establish the stationary turbulent motions) we have to 
consider turbulent fluctuations time-de-correlation which 
can be accounted for via introducing the characteristic 
function f(r](k),T) (with r]{k) the autocorrelation func¬ 
tion), [ 93 : 


0^2 _ 

f{r]{k),T) = exp[- (10) 

with To - the largest turbulent eddy turn-over time, K = 
fc/fcp. In this case, to determine the magnetic field two 
point correlation function in real (x) space, we have to 
account for magnetic field fluctuations at different time 
moments and at different positions, i.e. {bi{x.,t)bj{x + 
R, t -I- r)). Accordingly, in Fourier space, the two-point 
correlation function will be determined by the F^(k, t) 
with the time dependent spectral functions SM{k,t) and 
AM{k,t): 


^f(k,i) 


Pij{'k)SM{k,t) +mjikiAMik,t) 


X f{r]{k),t) 


( 11 ) 


Comparing with the stationary spectrum, see Eq. ®, 
we see that formally we replace = Fj^(k, t), by 

Fl^fkjt) = Fl^fk,t)f{r]{k),T), To avoid a complex de¬ 
scription of accounting time dependence of S{k,t) and 


A{k^t), we use the Proudman argument for kinetic tur¬ 
bulence ( 5 ^ , according which the description of decaying 
turbulence lasting for tt can be replaced by the descrip¬ 
tion of stationary turbulence with time duration o/tt'/ 2. 
Below we briefly discuss our approach. 

Turbulence during PTs generated through magnetic 
helicity can described through two major stages [s^ : 
during the first stage the main process is determined by 
the magnetic energy direct cascade that last few largest 
eddy turnover times tq = 27r/(fcoWo), where wq < 1 is the 
turbulent eddy velocity (vq ~ M for kinetic turbulence 
where M is the Mach number and vo — va for mag¬ 
netic turbulence) determined by PT and magnetogenesis 
model parameters, see Ref. [^, i.e. tt = sqTo (with 
So = 3 — 5). The magnetic field induces vorticity fluc¬ 
tuations, and at the end of the first (semi)equipartition 
between kinetic and magnetic energies is reached, that 
results in doubling the value of the source for GWs. The 
magnetic energy density power spectra are then deter¬ 
mined by the proper dissipation rate per unit enthalpy 
£m as: ^0 = with Ck constant order of unity, 

and £ = fco^Q. Note that the autocorrelation function 
rjfk) = [ 9 ^. Although the Kolmogoroff 

model is valid only for non-relativistic turbulence, while 
during PTs we might deal with uq — 1 (relativistic tur¬ 
bulence), our estimates for the amplitude and polariza¬ 
tion degrees of GWs signal are qualitatively justified, see 
[ 2 ^. The second stage consists on helicity transfer (in¬ 
verse cascade). The scaling laws for this stage are still 
under debate. Based on our previous consideration [s^, 
we assume that (i) the details of the scaling laws dur¬ 
ing this stage will not affect substantially our estimates; 
(ii) instead of considering decay-turbulence we will again 
consider the stationary turbulence with scale-dependent 
duration time. Then we obtain for the helical Kolmogo¬ 
roff model, Aq = Ckcc jfkocA 1 where a is the mag¬ 
netic helicity dissipation rate per unit enthalpy, leading 
to AojSo = ajieko). The helical Kolmogoroff model is 
mainly relevant for weakly helical fields, |A(fc)| <g; S'(fc), 
which is a case of magnetic fields generated during PTs.^ 
According to results of recent numerical simulations, see 
Ref. [ 5 ^, the weakly helical turbulence even account¬ 
ing for the free decay of turbulence, show an establish¬ 
ment of the spectra in a good agreement with the heli¬ 
cal Kolmogoroff model, as well as equipartition between 
magnetic and kinetic energy densities. Thus we adopt 
ns = —11/3 and nn = —14/3 for the inertial range, with 
So = and Aq = hSo with h which determines 

the fraction of helicity dissipation, h = a/{eko)- The 
turbulence fluctuation velocity vo = va is determined 
by the magnetogenesis mechanism and for the model of 
our interest is given by vq — 0.2 (BeSjn — 5 • 10“^ G) 


^ Note that the substantially helical case is usually determined by 
the helicity transfer (inverse cascade), So = CscPF and Aq = 
Cyin2/3 [g|. 
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[13 for EWPT model of Refs, [zl ^ and no — 0.01 
~ 2 • 10-® G) in for QCDPT model of Refs. 


IV. GRAVITATIONAL WAVES SIGNAL 
AMPLITUDE AND POLARIZATION 

In this section we compute the GW signal (stain ampli¬ 
tude) and polarization degree from the hydro and MHD 
turbulence generated during the first order cosmological 
PTs. To determine the amplitude of GWs we proceed as 
it is described in Ref. [13 ■ We derive the energy density 
spectrum of the GWs at the end of PT (in our approx¬ 
imation the end of turbulence). The energy density of 
GWs is given through the ensemble average as 

PGw{x,t) = (12) 

As we noted the rescaling of the GWs amplitude and 
frequency given through Eq. © is irrelevant when com¬ 
puting the polarization degree, while it is crucial for es¬ 
timation of GWs energy density. 


A. Gravitational Wave Signal 

Assuming the homogeneous and isotropic turbulent 
source lasting for tt, and using far field approximation, 
see [13|, the total energy density of GWs at a given spa¬ 
tial point and a given time can be obtained by integrating 
over all sources within a spherical shell centered at that 
observer, with a shell thickness corresponding to the du¬ 
ration of the turbulent source (in our case the duration of 
PT), and a radius equal to the proper distance along any 
light-like path from the observer to the source (causality 
requirement), and then 

PGw(w) = = 16Tr^uj^Gw'^TTH,jij{uj,uj), (13) 

a In O' 

where ui is the angular frequency measured at the mo¬ 
ment of generation of GWs, and is a compli¬ 

cated function of oj (which is computed through using of 
aero-acoustic approximation and Millionshchikov quasi¬ 
normality M) , given as 

- pkn JL 

“p(-;^)“f<-? 7 iF 73 ) (I'l) 

Here, erfc(a;) is the complementary error function defined 
as erfc(a;) = I — erf(x), where erf(a;) = dy exp{—y'^) 
is the error function (99| - The integral in Eq. dMD is 
dominated by the large scale (fc ~ fco) contribution so, 
for direct-cascade turbulence during the first stage (direct 
cascade, see Sec. Ill), the peak frequency is 

c, koM. (15) 


where M is Mach number. To compute the GWs sig¬ 
nal arising from the inverse cascade stage we have con¬ 
sider two models separately: Model A assumes that the 
correlation length during the inverse cascade scales as 
oc and Model B corresponds to the correlation 
length scaling as oc We obtain that in both 

models peak frequency during the second stage are equal 
and are determined by the Hubble frequency as [s^ : 

^ 2nH, (16) 


while the GW amplitude are slightly different in Model 
A and B as 


(k, uj) ~ H,jij (0, w) = 


exp 




r^°dk 

127r3/2fco Jks 


erfc 




(17) 


and 


iJy,,(k,w)(A) - i7y,,(0,w) = 


exp 


C*M2fc3 ) 


6^3/2fc3/2 fc7/2 


erfc 


cy^Mfc 3/2 


.(18) 


Here determines the amount of initial magnetic he- 
licity and is equal to C* = (a(x) • b(x))/(,^M^^M) (with 
a(x) = A(x)/w normalized vector potential), and ks = 
2 t:/I s is the typical scale at which the inverse cascade 
stops: either because the cascade time Teas reaches the 
expansion time scale or because the characteristic 
length scale Is — reaches the Hubble radius 
The value of ks can be found by using above conditions, 
being equal to ks = koC*^^'^i'y/My^'^. Note that the 
integrals in expressions Eqs. (ITT)) - (TTSl) are In this dom¬ 
inated by the large scale k ^ ks contributions and are 
maximal at 


B. Gravitational Wave Polarization 

To compute the polarization degree of GWs we need to 
estimate the tensor perturbations source two point corre¬ 
lation function r) = (n^^p*(k,t)n|y (k',t'-|-r)), 

which can be expressed through the forms and 

■A-ijim (which are defined below Eq. ([51)), as 

= (27r)3<5(3)(k_k/) 

X [M^jlmS{k,T) + iAijlmQ{k,T)] . (19) 

As we discussed above for the non-stationary turbulence 
S{k,T) and A{k,T) are complex functions of r and k 
(see Eqs. (flUl) - (fTTl) for the time de-correlation func¬ 
tion and the magnetic field spectrum). Following Ref. 
[^ we split the spatial and temporal dependence as 
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S{k,T) = S{k)Ds{T) and A(fc, r) = A{k)DA{T) which is 
a valid approximation for k k^ (the range which mostly 
contributes to the GWs signal). Generalizing the station¬ 
ary case @ through accounting for the time-dependent 
functions Ds{t) and Da{t) the forms for 5(fe,r) and 
Q{k,T) are given as, 


S{k,T) = -^^^Jd^piJd^p2S^^^k-pi-p2) 

X [{l + a^){l + /3^)Dl{T)S{pi)S{p2)+ 

+ Aa/3DAiT)A{pi)A{p2)], (20) 

w'^Ds{t)Da{t) 


Q{k,T) = 


d-^pi / d'^p2 


1287r6 

X - Pi - P2) [(1 + a^)j3S{pi)S{p2) 

+ {.^ + P^)aA{pi)A{p2)] , ( 21 ) 


where a = k • pi and (3 = k • P2 with pi = pi/pi 
(Pi = |Pi) and P2 = P2/P2 (pi = Ipi). The helical 
source term Q{k,T) vanishes for turbulence without he- 
licity. Since in the helical Kolmogoroff model the time 
de-correlation is mostly determined by the energy den¬ 
sity dissipation, thus at first order approximation we can 
assume that Ds{t) ~ Da{t) (both functions are mono- 
tonically decreasing functions). By next we should con¬ 
nect S and Q with the H{k,t) and 'H(fe,t) functions, 
which determine the GWs polarization degree, see. for 
details Ref. [32|. 

Magnetic helicity generated via bubble wall collisions 
during the first order PTs is determined by the corre¬ 
sponding energy scales Apt (Apt — 100 GeV and 0.15 
GeV for EWPT and QGDPT respectively); and the PT 
bubble lengths (/b). In addition, the bubble wall velocity 
substantially affects the turbulent motions development 


10 sec) is equal to va — 0.2, and assuming around 
100 bubble within Hubble length scale, fractional helicity 
is — 0.01. Gomoving magnetic helicity itself is 

expressed: 


^(EW) ^ 

n-MA - i25GeV ’ 


(23) 


Assuming that the magnetic field (with va — 0.01) is 
correlated over the wall thickness (the QGD momentum 
0.15 Gev) [o^, results in extremely small magnetic he¬ 
licity generated during the first order QGDPT (at time 
moment ~ 10“® sec). 


a/(QCD) 

n-M 


(G) 


('^(QCD)^2 

0.15GeV 


(24) 


which corresponds 1. On the other hand, mak¬ 

ing the field correlated over the bubble length scale, leads 
to the fractional helicity ~ 0.2 

Using the approximation given above, the polarization 
degree of GWs, {k), for the Kolmogoroff helical tur¬ 
bulence model can be estimated through (in our simpli¬ 
fied description the time dependence is canceled because 
of Ds{t) ~ Da{t)): 


-pGW^k) 


n{k) 

H{k) 


Ta{K) 

MK) 


(25) 


where K = k/kg is a normalized wavenumber, and 

Is{K) ~ JdPi PiJdP 2 P20 [(1 + a^)(l + Pl)Pi^P 2 ^ 
+4h^apPpP^^P^^] , (26) 

Ia{K) ~ 2h jdPi Pi jdP2 P 2 & [(1 + aDPpP^^P^^ 


In the present paper we focus on magnetic helicity 
generation mechanisms following Refs. In the 

framework of these magnetogenesis scenarios magnetic 
helicity during PTs with the magnetic wall in the x — y 
plane is given by 

Um = A,B, or Um = , (22) 

Apt 

We note that the EWPT energy (mass) scale is approx¬ 
imately equal to Higgs mass, Aewpt — Afniggs- The 
model parameters such as bubble and wall sizes and wall 
velocity depend on PT modeling, and are determined in 
Refs. uM- We quote also the physical values of 

magnetic field amplitudes as ~ 6.45 x lO^GeV^ 

[iilll and ~ 1.5 X IQ-^GeV^ [1^. 

The fraction of initial magnetic helicity (G) can be 
expressed in terms of the magnetic field correlation length 
(which can be taken to be equal to 4) and the maximal 
allowed leiigth scale as G — Following 

Refs. uM normalized magnetic field generated 
during the first order EWPT (at the time moment ~ 


+(l + /?2)apPi"-p"1 • (27) 

Here h is the model parameter (which is related to the 
helicity fraction, see below, and we assume it to be equal 
to 1, 0.5, 0.1. Pi = pi/fco, P2 = P2/fco, oip = (/G -I- 
P2 - P|)/(2KPi), /3p = (P^ + P| - PI)H2KP2), 0 = 
6»(Pi-bP2-A:)6»(Pi-bA:-P2)6»(P2-f AT-Pi), and 6» is the 
Heaviside step function which is zero (unity) for negative 
(positive) argument. The integration limits ranges from 
1 (we discard the source existence for the wavenumbers 
below fco) to ku/kQ. 

We emphasize that the fractional helicity parameter 
G discussed above is defined through normalized mag¬ 
netic helicity (the integral quantity), while the parameter 
h = a/{koe) is defined through the normalized magnetic 
energy density and normalized magnetic helicity (e.g. it 
is determined by the power spectra for magnetic energy 
density and helicity at small length scales). Under the 
model adopted here (the Kolmogoroff helical turbulence 
with ns = —11/3 and ha = —14//3) these both quanti¬ 
ties coincide C,* — h. 

To keep our description as general as possible we 
present our results for the GWs polarization degree 
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FIG. 1: The GW polarization degree {K) {K — k/ko) in terms of the model parameter h — 1.0, 0.5, 0.1 


{k) in terms of the normalized wavenumber K. As 
we discussed above the typical wavenumber ko is deter¬ 
mined by the turbulent eddy length scale Iq (fco = 27 r//o) 
and is significantly different for EWPT and QCDPT. The 
model parameter h, which determines the helicity frac¬ 
tion, varies depending on the magnetogenesis model. The 
results for {k) with h = 0.1, 0.5, and 1 are shown in 
Figure 1. 

V. CONCLUSIONS 

We computed the GWs signal produced during first 
order cosmological PTs through hydro and MHD turbu¬ 
lence. We also derive the polarization degree of GWs 
assuming the validity of the helical Kolmogoroff model, 
shown in Figure 1. The GWs polarization is present at 
the background level, and for maximally helical sources 
the polarization degree approaches unity at its maxi¬ 
mum, around k ^ 2fco, and decreases fast at small scales 
k ^ ko- The formalism presented in this paper might be 
used to estimate the polarization degree of GWs from he¬ 
lical hydro and M HD t urbulence in the differ entia l rotat¬ 
ing neutron stars [io3 or stellar convection |l0l| . Note 
from Figure 1 that the detectability of the polarization 
degree is determined by the helicity fraction parameter 
h. We plan in our future research to make estimates 
of h values depending on magnetogenesis models during 
cosmological PTs. 

Previously we have estimated the GWs amplitude, 
hc{f), from the first order FWPT and QCDPT [sS, 
thr oug h assumptions of non-helical magnetic fields [60, 
III, [^. We have shown that EWPT generated GWs are 


potentially detectable through LISA-like missions in 
the case for strong enough EWPT (for QCDPT gen¬ 
erated GWs detection prospects see ReF. [iO)- In the 
present paper we expand our previous results by consid¬ 
ering GWs from helical magnetic and hydro turbulence. 
Probing the circular polarization of GWs background is 
a challenging task [^, and it is quite difficult at the 
monopole mode. To detect the circular polarization at 
the dipole or/and the octopule mode requires at least 
the system of two unaligned detectors, and LISA was 
designed ideally to provide detection of these anisotropic 
components whose magnitudes are small as 1% of the de¬ 
tector noise [6ll-l^. The similar detection prospects are 
expected from eLISA data. The planned eLISA mission 
[491 is originally design to detect un-polarized GWs bank- 
grounds (including GWs from cosmological PTs) |l02j | 
with the sensitivity given on Fig. 1 of Ref. [5l| (see also 
comparison with EISA’s sensitivity). Although the GWs 
polarization detection is beyond the currently discussed 
eLISA science, our study should help further develop¬ 
ments. 
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